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Abstract. We study the bipartite entanglement of two general classes of 
heterogeneous spin-(l, ^) and homogeneous spin-1 systems. By employing the spin 
correlation functions, we obtain the reduced two-spin density matrix (DM) and the 
£h ■ negativity for these two classes of quantum spin models. We show explicitly that 

in addition to the one and two-point correlations, the triad and quad correlations 
(t*^ = (S a Spss) and q^L = (SaSpsss^) where a, /3, S, 7 = ±, z) play crucial role in the 
bipartite entanglement between spins s > 1/2. These correlations represent the spin i- 
quadrupole and quadrupole-quadrupole correlations, respectively. These correlations 
' do not appear in the spin-i models. Our results are general and applicable to the 

. different several models of interest with higher reflectional, translational, spin-flip and 

U(l) symmetries. The entanglement of many attractive models are investigated. 
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1. Introduction 

In many-body physics, quantum spin systems have received a lot of attentions due to 
their importance in the description of different phenomena such as antiferromagnetic 
Mott insulating behavior, Bose-Einstein condensation (BEC) [1], spin supersolidity [2] 
and gapless Liittinger liquid physics. [5] Apart from pure scientific interest, spin systems 
possess a potential for various technological applications such as nano-technology and 
quantum computations. They are prototype realization of many relevant properties of 
quantum implementation in advance technologies such as fabrication of spin transverse 
torque switching devices, jl] Therefore different aspects of a quantum phase are of utmost 
importance for scientists and engineers. 

Entanglement is a kind of quantum nonlocal correlation which is used as an 
indicator of quantum phase transitions. [HI E] The collective behavior of a system 
is associated with the development of this nonlocal correlations as well as classical 
correlations. Moreover condensed matter systems can be efficiently simulated using 
techniques related to entanglement. [7] The usefulness of entanglement in condensed 
matter physics leads us to study it in the parity invariant (PI) quantum spin systems. 
Recently, the geometrical entanglement has been studied for one- dimensional systems 
with space or space-time inversion symmetries. [8] The main reasons to consider the 
systems with parity symmetry are as follows: (i) The external parameters such as single 
site and lattice anisotropics, spin-orbit couplings and external magnetic and electric 
fields do not break the parity and we may investigate the simultaneous effects of these 
parameters in the entanglement properties of the state, (ii) It is a general model and the 
results are employed to study the entanglement properties of many interesting systems 
with higher reflectional, translational, spin-flip and U(l) symmetries. 

In the present paper, we investigate the bipartite entanglement in parity 
invariant (PI) heterogeneous mixed spin-(l, |) and homogeneous spin-1 systems. Many 
theoretical and experimental investigations of the physical properties of S > \ spin 
systems have been done since Haldane's conjecture[9] and synthesis of bimetallic 
molecular magnets and supramolecular composite systems. [10J An interesting system 
with parity symmetry is the XXZ spin-(a, |) or spin-1 model in a transverse magnetic 
field with Dzyaloshinskii-Moriya interaction arising from spin-orbit 

coupling). We obtain an expression for the two-spin reduced density matrix (DM) 
and the negativity in terms of the spin correlation functions. All the information needed 
to analyze bipartite entanglement is contained in the two-spin reduced density matrix, 
obtained from the wave function of the state after all the spins except those at positions 
i,j have been traced out. The resulting reduced density matrix represents a mixed 
state of a bipartite system for which a good deal of work has been devoted to quantify 
its entanglement. The spin correlation functions show the collective critical behaviors 
of particles, as well as their ordering and fluctuations in a system. Utilizing spin-spin 
correlation functions, enables us to investigate the effects of various external parameters 
(magnetic and electric field, anisotropies and temperature) to study the entanglement 
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properties of a state such as the dynamics of entanglement, entanglement witnesses 
and quantum fidelity. We show explicitly that in addition to the magnetizations (S?) 
and two-point correlations (S^S"), the information of triad and quad correlations 
(ta/3 = (SaSpSg) and q s ^p = (S a S/3Sss^) where a, ft, 8, 7 = ±, z) play key role in the 
bipartite entanglement in ferrimagnetic (S = l,s = |) and antiferromagnetic S = 1 
spin systems. They are indeed the spin |-quadrupole and quadrupole-quadrupole 
correlations. Finally, employing our results we study the entanglement in different 
several models of interest with higher symmetries. 

2. PI heterogeneous (S = l,s = |) spin systems 

A parity invariant state satisfies the relation [V, p] — where the parity operator V is 
defined as: 

V = c -*rE«W+*7> = - (<?f) 2 )4 (1) 

For the homogeneus spin-| system the parity operator is reduced to the more familiar 
form: Yli s f • Let us utilize spin correlation functions to obtain the elements of DM. A 
relationship between concurrence [UJ and two-point correlation functions of the spin-| 
anisotropic Heisenberg XY model in the presence of a transverse magnetic field has 
been already presented. [12] For any spin-| system, all elements of DM are acquired 
by using the magnetizations and the two-point correlations. However, these functions 
are not sufficient to attain the DM of PI spin-(l, |) or spin-1 systems. As a matter of 
fact the quadrupolar effects of the spin-1 particles play important role on the physical 
behavior of the system and we should employ them to achieve the DM. Defining the 
triad correlations as t s aj3 = (SaSpSs) where a, (3,8 = ±, z and working in the basis of 
factorized state the DM of a PI state is given by 
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where 

4 = (S a S p [l ± 2s z }) = (S a Sp) ± 2t z af} , 

t& z ±p = ([! ± S Z ]S P 8 S ) = g p& ± 1? zl3 , gp S = (Sps 5 ). 

Some off-diagonal elements in (j2J) like t~± = (S Z S±S-) and t± z = (S±S z s + ) are the 
correlations of the transverse component of spin-i with operators (S Z S±) and (S±S Z ). 
These are nothing but the quadrupolar operators [13] which have been introduced to 
describe the quadrupolar order in the spin-1 bilinear-biquadratic Heisenberg model. [H] 
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The reduced density matrix strongly depends on the correlations between the spin-| 
particle and the quadrupoles. Some of these correlations are suppressed at the presence 
of other symmetries. For example, for states with parity and spin- flip symmetries |15] the 
equalities 3lt+ z = 3ltt z = (3? denotes the real part), t+ + = t z __ and M z = m z = t z zz = 
are gained. (M z = (S z ) and m z = (s z ) are the sub-lattice magnetizations in z direction). 
Moreover, adding U(l) symmetry to the parity, the equalities t z __ = t z __ = t z _ = 
are obtained. By making use of vigorous analytical and numerical approaches such as 
bosonization, Green's functions techniquw, exact diagonalization Lanczos and density 
matrix renormalization group methods we can study the behaviors of the correlation 
functions. 

Negativity for spin- (1,1/2) system 
Utilizing the reduced density matrix (j2j) we can obtain the entanglement, as quantified 
by negativity and relative entropy of entanglement, between two spins in a PI spin-(l, |) 
state. In this paper we focus on the negativity and the relative entropy of entanglement 
is left for the future. Negativity is based on the trace norm of the partial transpose of 
density matrix and is given by the formula: 

AW) = llpTl ^~\ (3) 

which corresponds to the absolute value of the sum of the negative eigenvalues of partial 
transpose density matrix, p Tl .[T6] 

The negative eigenvalues are obtained by solving two cubic equations: x 3 + a±x 2 + 
b±x + c± = where a± = —\ ± (m z — 2t z zz ) are bounded by — 1 < a± < and 

o± = ^£A(M, ± 2g zz f - [tf z f + |*it| 2 ] 

+ ± + t:+(t z ±tt + =f t z z U+_)i 

b± = ^{tfAKt + < +z -) ~ (M z ± 2g zz f - \t z _y - 8[|t; + | 2 + |t;_| 2 ]}. 
By defining the parameters f± = |a| — ^b±, r± = \{b±a± — 3c±) — and 

_3 

cos^i = r±f± 2 the standard solutions of the cubic equations are 

± r, l~T~ ___I®±\ a ± 



x i = 2V/±cos v , 

± rr ,9± + 2vr a± 
x 2 = 2 V/± cos( ) - 

4 = 2v / Sco 8 (^) - f , (4) 

where 6± G [0,7r]. The negativity is the absolute value of the sum of the negative 
x's. All possibilities for the sign of eigenvalues have been explicitly shown in Fig. (JTJ). 
If (xf,x^) G region (1), then A/12 = ^1=2 z\ x t \ 2 + \ x i\ 2 an d the state is maximally 
entangled. However in region (3) negativity is exactly zero which corresponds to a fully 
factorized state and the state is partially entangled in other cases. A system with parity 
symmetry is the ferrimagnetic XXZ spin-(l, |) model in a transverse magnetic field with 
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Figure 1. Different possibilities for the sign of the solutions of the two cubic equations. 
By Comparing —3-^/f± and — 6^//± with 1 we find the allowable regions simply. 



Dzyaloshinskii-Moriya interaction The parameters 9± and f± depend 

on the transverse magnetic field and anisotropy parameter D. Quantum phases are 
specified by D and magnetic field and the negativity would give useful information on 
the entanglement of different phases of the system. 

To see more advantages of these results it is worth to investigate the entanglement 
of states with higher symmetries. By adding the spin-flip symmetry to the parity, the 
coefficients of cubic equations are simplified as a = a + = a_ , b = b + = , c = c + = c_ 
and the negative eigenvalues have been indicated in Fig.fJTJ. For b < 0, the state 
is always entangled and the entanglement depends on 9. For small values of 9 the 
state could be maximally entangled and the entanglement is decreased by increasing 
9. At b = the entanglement is independent of 9 and the state is always partially 
entangled. For < b < -k and small values of 9 the state is fully factorized however 
it is partially entangled by increasing 9. Finally for b > A the state is always fully 
factorized. Some examples of models with such symmetries are as follows: (i) Mixed 
spin XXZ and XY systems in a transverse magnetic field [17J, and (ii) Mixed spin Ising 
model with (or without) next nearest neighbor interaction in the presence of a transverse 
and a longitudinal field. Ferrimagnetic Ising models have been studied with a renewed 
interest, especially in connection with compensation points. These points occur at 
temperatures below the critical temperature, at which the sublattice magnetizations 
cancel exactly, and the total magnetization changes sign. These systems maybe used in 
technical applications, most notably in magnetic recording. 

Finally a mixed state with spin-flip, £7(1) and parity symmetries is entangled if 
|[1 -\{Sf) - g zz - \J[\{S 2 Z ) + g zz - l] 2 - 8[t; + ] 2 ] < 0. All anisotropic bond-alternating 
spin-(l, |) ferrimagnets (with dimerization, trimerization and tetramerization) in a 
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longitudinal field, the ferrimagnetic XXZ models in a longitudinal field have these 
types of symmetries. The results are applicable to the synthesized ferrimagnetic 
compounds with Ising-like anisotropy such as C oNi(E DT A)6H20 where (S, s) = 

(i(Ni),i(co))m 



3. PI homogeneous S — 1 spin systems 



In this section we obtain an expression for the reduced density matrix of PI spin- 
1 systems. The parity operator is simplified as V = Ylii^ ~ 2(<S!f) 2 ). In addition 
to the sublattice magnetizations and two-point correlations other correlations should 
be employed to find all elements of the DM. By introducing the quad correlations 
1% = (SaSpSss^) (indices are ± and z) and working in the basis of factorized state 
the DM is given by 
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where \S\ — \s\ — 1 and 

gS ± = (5 Q ^(l±^)) = C±gS, 

g^ 2± = (5 a (l ± 5 Z ) S5S7 ) = # ± g£. (6) 

Appearance of the quad correlations like gf- = (S z S-S z s±) and g^f = (S z S-S±s x ) 
in Eq.([5]) is noticeable. By writing S^SL and s z s± or s±s 2 in terms of quadrupolar 
operators [13] we find that the quadrupole-quadrupole correlations play an essential role 
in the DM of homogeneous spin-1 models. This kind of correlations appear neither 
in heterogenous spin-(l, |) systems nor in homogeneous spin-| models. They are the 
characteristic behavior of the spin-1 systems. Adding further symmetries to the parity, 
suppresses some of such correlations. 

Comparing the DM of three spin-|, spin-(l, |) and spin-1 systems, we conclude 
that always one component of spin-| operator, and two components of spin-1 operator 
contribute to the triad and quad correlation functions. Moreover, all information for 
the bipartite entanglement of a spin-| system are gathered from the magnetizations, 
staggered magnetizations and the two-point correlations, however in a heterogenous 
spin-(l, ~) and a homogeneous spin-1 system in addition to one and two-point 
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correlations the triad and quad correlations play key roles in the bipartite entanglement 
of PI states. 

Negativity for spin-1 system- For systems with dimension higher than 2 <g> 3, 
inseparable states with positive partial transpose (PPT) exist, i.e. the PPT or Peres- 
Horodecki criterion is in general a necessary but not a sufficient condition. [19] Therefore 
a state is certainly entangled, whenever one of the eigenvalues of p Tl is negative. 

By solving a cubic and three quadratic equations the eigenvalues of p Tl are obtained 
in terms of magnetizations, two-point, triad and quad correlation functions. The 
expressions are too lengthy and we have not brought them here. However their forms 
are simplified by considering other symmetries. 

States with translation, reflection and spin-flip symmetries - At the simultaneous 
presence of the translation, reflection and spin-flip symmetries only the following 
eigenvalues could be negative. 

xi = K + ±G„, x 2 = K_±G + , x 3 = -(r 1 / 3 + ^-), 
1 



x 4 = - (^M_ + P_ - y (M_ - P_) 2 - 16(3g-+) 2 J . (7) 

where r'(< 0), r, N±, K±, G± and P± are functions of quad correlations (see appendix) 
and S denotes the imaginary part. The anisotropic XXZ spin-1 models in a transverse 
magnetic field have these symmetries. By varying the non-commuting magnetic field 
the system experiences different phases in which the correlation functions and thus the 
negativity behaves differently in each phase. 

States with parity and £7(1) symmetries - For states with parity and £7(1) 
symmetries, among 6 eigenvalues of the three quadratic equations only the following 
eigenvalues 

1 



1 

Xl = 4 



A^-JAi + 4\ql 



x 2 = - 



B_ - J Bl + 4|g£l' 



could be negative, where A± = [q zz + + ± (f z - z +\i B± = ^ ( f Z z- + \- ^ ne °^her 

eigenvalues are given by the standard solution of the cubic equation with coefficients a, 
b, and c (see appendix) where a is bounded as — 1 < a < and the conditions to have 
a negative eigenvalue have been discussed in Fig.([T]). 

States with parity, spin-flip and U(l) symmetries - Finally, for PI states with 
additional spin-flip, 27" (1) and translational symmetries, the following eigenvalues could 
be negative: 

x l = l -( K D + ±s]Dl + ?,2\qtl}^ X2 = ^~, x 3 = K ± . (9) 

where D± = P + ± ^q z z X z z ~ ■ These results give the negativity of the anisotropic XXZ 
spin-1 models in a longitudinal magnetic field. Adding a single site uniaxial anisotropy 
to this model does not break the symmetry and the negativity is given by the negative 
eigenvalues of ()9]). The anisotropic XXZ spin-1 in a longitudinal magnetic field with 
single site anisotropy is used to describe the spin super-solid and spin-flop phases. [2J [20] 
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Meanwhile NiCl^ — 4SC ( N H<>)<)( DT N)^T\ is a quantum spin-1 chain with strong easy- 
plane anisotropy and a new candidate for the BEC of the spin degrees of freedom. 
Finding the negative eigenvalues of the entanglement of the BEC state is determined. 



4. Summary 

In this paper we have studied the bipartite entanglement, as quantified by negativity, 
between two spins in two general classes of spin system, the parity-invariant (PI) 
heterogeneous spin-(l, |) system and parity invariant homogeneous spin-1 models. All 
the information needed to analyze bipartite entanglement is contained in the two-qubit 
reduced density matrix, obtained from the wave function of the state after all the spins 
except those at positions i,j have been traced out. The resulting reduced density 
matrix represents a mixed state of a bipartite system for which a good deal of work has 
been devoted to quantify its entanglement. To see the effects of collective behaviors of 
particles on the bipartite entanglement we have expressed the reduced density matrix 
and the negativity in terms of spin correlation functions. We have shown explicitly that 
the triad and quad correlations like (S z S±s z ) and (S z Ss z s±), which demonstrate the 
quadrupole-single spin and the quadrupole-quadrupole correlations, are emerged in the 
reduced density matrix and the negativity. We have also found that in mixed spin- 
(1, |) and spin-1 systems, in addition to magnetizations and two-point correlations, the 
information coming from the quadrupolar operators (S Z S± and S±S Z ) play an essential 
role in the pairwise entanglement. These correlations do not appear in the spin-| models. 
Utilizing our results it has been also achieved the entanglement of several models of 
interest with higher reflectional, translational, spin-flip and U(l) symmetries. 

The effects of the triad and quad correlations on the other useful quantities such 
as entanglement entropy, quantum fidelity will be investigated in future. 
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6. Appendix 



N ± = qtt±g zzi K ± = ~[qf/ + ±q*+], G ± =i[£r±(£ 
P_ = N + - qtt, M_ = AL - qZZ, r' = ~(L+ + Aq " 



Z T z \ 
z+z~ ) ' 



1 

96 



1 



-{P + + 4qf +z Z)M + - qf + *ZP + + 2(fe£] a + [q~*] 2 ) 



4k 



z + z~]2 

z+z-l 



I 2 

zz+\ 



Entanglement in heterogeneous spin-(l, \) and homogeneous spin-1 systems 



9 



Cubic equation: 

ax 3 + bx 2 + cx = 0, 



b = ^(C+C- - kill 2 + - ltf-l 2 D. 
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